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ABSTRACT 


The  relation  between  natural  electrical 
currents  flowing  in  the  earth,  or  telluric 
currents,  and  subsurface  structure  was  studied 
using  procedures  developed  by  previous  investi¬ 
gators.  The  scale  of  the  study  was  larger  than 
had  been  attempted  previously,  and  the  results 
were  not  as  conclusive.  It  was  found  that  the 
telluric  currents  are  limited  in  their  application 
to  subsurface  exploration  by  the  non-unif ormity  of 
the  telluric  field  over  very  large  areas.  This 
practical  limit  of  telluric  exploration  is  thought 
to  be  several  hundred  miles. 

A  second  method  of  subsurface  exploration 
was  studied  using  electromagnetic  induction  rel¬ 
ationships  between  telluric  currents  and  the 
earth* s  magnetic  field  at  the  surface  of  the  eaasth. 
The  relation  between  electromagnetic  induction  and 
subsurface  structure  had  been  predicted  theoretic¬ 
ally  but  no  experimental  results  have  been  pub¬ 
lished.  It  was  found  that  measurements  of  the 
necessary  order  of  accuracy  were  possible,  but  that 
the  method  cannot  be  thoroughly  tested  until  new 
procedures  are  developed  for  analysing  the  data.. 
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CHAPTER  ONE 


INTRODUCTION 

The  natural  electromagnetic  field  of  the  earth 
is  affected  at  the  surface  by  the  electrical  pro¬ 
perties  underground.  It  is  the  purpose  of  this 
research  to  investigate  experimentally  two 
particular  methods  of  subsurface  exploration  using 
the  earth* s  electromagnetic  field. 

The  investigations  of  each  of  the  two  methods 
of  exploration  are  developed  in  a  parallel  manner 
throughout  the  thesis.  The  first  method  investi¬ 
gated  is  that  of  Telluric  Prospecting,  This 
method  uses  the  variations  in  the  electric  currents 
flowing  in  the  earth  —  "telluric  currents"  —  to 
determine  subsurface  structure*  The  second  method 
is  called  "magneto-telluric  prospecting"  and  uses 
the  relation  between  the  electric  and  magnetic 
components  of  the  electromagnetic  field  at  the 
surface  for  determining  subsurface  structure. 

Telluric  prospecting  has  been  attempted  in 
Europe1*2  and  the  United  States. The  scale  of 

i 

Schlumberger,  Trans.  Am.  Geoph.  Union.  3:271 

(1939) 


the  prospecting  was  comparatively  small  —  less 
than  one  hundred  miles  across  the  area  covered. 

The  results  obtained  showed  a  correlation  between 
the  telluric  currents  at  the  surface  and  subsurface 
structure.  The  present  application  of  the  method 
of  telluric  prospecting  was  applied  over  a  larger 
scale  (about  one  thousand  miles)  and  was  an  attempt 
to  correlate  telluric  currents  across  the  Canadian 
Prairies  with  the  depth  of  the  Precambrian  basement. 
Very  few  details  have  been  published  in 
English  about  methods  of  telluric  prospecting. 

Those  which  are  available  are  concerned  with  local 
structural  features.  It  is  the  purpose  of  this 
study  to  give  details  about  techniques  in  telluric 
prospecting  and  to  apply  these  techniques  to  a  large 
scale  structural  feature.  It  was  further  hoped  that 
a  large  scale  study  might  yield  more  information 
about  the  nature  of  telluric  currents  in  general 
over  a  large  area. 


Schlumberger  and  Kunetz,  C .  B.  Acad.  Scl, 
Paris.  223:551-53  (1946). 
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Dahlberg,  Geophysics .  10:494  (1945), 

4 

Boissannas  and  Leonard on,  Geophysics . 
13:38?  (1948) 
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No  results  of  magneto-telluric  prospecting 
have  been  published  up  to  the  present.  Therefore, 
it  seems  desirable  that  an  experimental  investigation 
of  magneto-telluric  prospecting  be  attempted  follo¬ 
wing  the  outline  of  the  theory  developed.  The  theory 
has  been  developed  by  Cagniard,^  who  suggested  that 
subsurface  structure  could  be  determined  from  phase 
and  amplitude  relationships  between  the  components 
of  the  telluric  and  magnetic  fields  at  the  surface. 

The  depth  of  the  Precambrian  basement  was 
used  as  a  known  geological  feature  in  the  magneto- 
telluric  prospecting  as  in  the  telluric  prospecting. 
Experimental  stations  were  established  at  locations 
which  gave  a  wide  range  of  depth  of  the  Precambrian 
basement.  The  results  of  the  magneto-telluric 
exploration  were  compared  with  the  known  geologic 
structure  to  determine  the  success  of  the  method 
as  applied. 

In  the  remainder  of  the  thesis,  the  corres¬ 
ponding  features  of  each  method  are  outlined  in  a 
parallel  manner.  The  following  three  chapters 
describe  the  theory  of  the  electromagnetic  field, 
the  details  of  the  field  methods  used,  and  the 
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Cagniard,  Geoichvsics .  18:605,  (1953) 
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methods  used  in  analyzing  the  data  obtained.  Typical 
cases  of  each  method  of  prospecting  are  analyzed  in 
Chapter  Five,  and  the  results  of  the  overall  explo¬ 
rations  are  interpreted. 


CHAPTER  TWO 


THEORY 

a)  The  Earth* s  Electromagnetic  Field 

The  magnetic  field  of  the  earth  is  constantly 
varying.  It  can  be  broken  up  into  two  portions: 
the  ambient,  or  dipole  magnetic  field,  which  varies 
only  slowly  over  the  course  of  many  years,  and  the 
varying  magnetic  field,  composed  of  both  irregular 
random  variations  and  a  number  of  periodic  regular 
variations , 

The  "earth* s  electromagnetic  field"  as  used 
here  refers  to  the  varying  magnetic  field  plus  the 
varying  electric  currents  which  flow  in  the  earth. 
The  electric  and  magnetic  fields,  considered  to¬ 
gether,  make  up  the  earth* s  electromagnetic  field. 

The  connection  between  the  telluric  currents 
and  the  magnetic  variations  becomes  evident  when 
the  source  of  the  electromagnetic  field  is  con¬ 
sidered.  The  location  of  the  source  of  such  fields 
may  be  determined  by  means  of  a  spherical  harmonic 
analysis.  Such  an  analysis  was  used  by  Gauss  to 


Chapman  and  Bartels,  Geomagnetism ,  p.606, 
Oxford,  19^0. 
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show  that  the  dipole  magnetic  field  is  of  internal 

origin,  A  similar  analysis  applied  to  the  varying 

magnetic  field  shows  that  it  is  about  two-thirds 

due  to  sources  which  are  located  outside  the  earth 

itself,  and  about  one-third  due  to  sources  inside 

7 

the  earth.  The  internal  portion  of  the  magnetic 
variation  can  reasonably  be  explained  as  due  to 
electric  currents  induced  within  the  earth  by  the 
external  portion  of  the  magnetic  variation.  The 
external  sources  of  the  magnetic  variation  are 
believed  to  be  vast  current  sheets  located  in  the 
ionosphere  of  the  earth. 

The  source  of  the  electromagnetic  field, 
therefore,  is  situated  at  a  great  distance  away  from 

o 

the  earth,  Cagniard  points  out  that  the  length  of 
the  wave  of  a  disturbance  is  enormous*  For  a  distur¬ 
bance  of  thirty  second  period,  the  wave  length  would 
be  nine  million  kilometers,  which  is  fourteen  hundred 
times  the  radius  of  the  earth.  Such  a  wave  would 
extend  over  the  whole  earth.  This  widespread  nature 
of  the  electromagnetic  disturbances  ensures  a 

7  Ibid . .  p.  684 

p 

Cagniard,  Handbuch  Per  Physik.  p,  46 2, 
Berlin,  1956 
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a 

uniformity  over^large  area  of  the  earth. 

The  varying  electromagnetic  field  may  be 
observed  on  the  surface  of  the  earth.  It  is  seen 
as  a  continually  varying  phenomenon  whose  variations 
have  a  wide  frequency  range,  and  which  persist  for 
varying  lengths  of  time.  The  direction  of  the 
varying  field  is  continually  changing.  Because  it 
is  expensive  and  complicated  to  measure  all  of  the 
independent  components  of  the  field  simultaneously, 
it  is  usual  to  measure  only  a  few  of  the  components. 
Therefore,  only  the  horizontal  components  of  the 
electric  and  magnetic  fields  at  the  surface  will  be 
considered  in  the  following. 

Two  types  of  prospecting  using  the  earth* s 
electromagnetic  field  will  be  considered  in  detail. 
The  first,  known  as  telluric  prospecting,  uses  point- 
to-point  variations  in  the  telluric  currents  to  give 
subsurface  information.  The  second  type,  known  as 
magneto-telluric  prospecting,  measures  the  magnetic 
and  telluric  components  simultaneously  at  one 
station.  It  is  thought  that  subsurface  information 
can  be  obtained  from  amplitude  and  phase  relation¬ 
ships  between  the  electric  and  magnetic  components. 
The  theory  and  application  of  each  of  these  methods 
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* 
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will  be  treated  in  turn 
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b)  Telluric  Current  Theory 

It  is  found  that  a  varying  voltage  exists 
between  two  electrodes  inserted  in  the  ground* 

This  varying  voltage  is  associated  with  sheets  of 
electric  currents  in  the  earth,  known  as  telluric 
currents.  It  is  usual  to  refer  to  the  measurement 
of  these  voltages  as  measurements  of  the  telluric 
current.  This  practice  will  be  followed. 

Telluric  currents,  when  simultaneously 
measured  over  a  part  of  the  earth's  surface,  show 
similarities  at  different  locations.  However, 
though  similar,  earth  currents  do  show  variations 
with  location  which  are  thought  to  be  due  to  sub¬ 
surface  electrical  properties.  If  this  is  true, 
then  simultaneous  earth  current  measurements  over 
an  area  should  reveal  indications  of  subsurface 
structure. 

To  measure  earth  currents,  two  conducting 
electrodes  are  placed  in  the  ground  at  distances 
apart  which  may  vary  from  a  few  hundred  feet  to 
several  miles.  A  voltage  recording  device  is 
placed  between  the  electrodes.  This  gives  a  record 
of  the  potential  difference  between  the  two  elect¬ 
rodes,  which  may  be  interpreted  as  a  record  of 
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' 
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voltage  gradient  in  the  direction  of  the  line 
joining  the  two  electrodes.  Voltage  gradient  is  a 
vector  quantity  with  both  magnitude  and  direction. 

To  determine  the  potential  gradient  vector,  it  is 
necessary  to  measure  two  components,  preferably  at 
right  angles. 

Telluric  currents  flow  in  large  sheets. 

These  sheets  are  uniform  over  a  large  area  at  a 
given  time,  but  are  constantly  varying  with  time, 
both  in  amplitude  and  direction.  Because  it  would 
be  far  too  expensive  to  measure  the  telluric  currents 
simultaneously  at  a  great  number  of  stations,  it  is 
usual  to  measure  the  telluric  currents  only  at  a 
fixed  base  station  and  at  a  field  station  which  can 
be  mobile.  If  all  the  field  station  records  are 
compared  with  those  at  the  base  station,  results 
can  be  obtained  as  if  the  telluric  currents  were 
measured  simultaneously  at  a  large  network  of 
stations.  The  electrical  properties  of  the  sub¬ 
surface  can  then  be  determined  for  the  area  covered 
by  the  survey. 

Assuming  that  the  earth  current  sheet  is 
continuous,  there  will  be  a  linear  relation  between 
the  telluric  vectors  at  the  base  station  and  the 


11 


field  station*^  Letting  the  components  of  the 
telluric  field  at  the  base  station  be  x  and  y,  and 
at  the  field  station  X  and  Y,  then  the  linear 
relation  may  be  written: 


X  =  ax  *  by 
Y  *  ex  +  dy 


(1) 


Equations  (1)  are  a  linear  vector  transf ormation 
a  *b 


The  determinant,  A  =  I  , 

led 


,  is  invariant  under  the 
selection  of  co-ordinate  systems*  That  is,  whatever 
the  orientation  of  the  axes  used  to  measure  the 
telluric  currents  at  a  station,  the  determinant  A 
will  be  the  same*  Thus,  A  is  associated  with  a 
particular  location  and  is  invariant* 

The  quantities  expressed  in  equation  (1), 

(viz.,  X,  Y,  x,  and  y)  represent  potential  gradients 
measured  between  the  electrodes*  Measuring  the  value 
of  the  telluric  field  accurately  is  difficult, 
however,  due  to  potentials  developed  at  the  electrodes 
These  “contact  potentials are  super-imposed  upon 
the  natural  telluric  potentials.  Since  it  is 
difficult  to  measure  the  value  of  the  contact 
potentials,  it  is  not  practical  to  try  to  measure 
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Porstendorfer ,  Tellurik.  p*51,  Berlin,  1954. 
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the  total  telluric  field  accurately. 

Instead,  measurements  are  confined  to  vari¬ 
ations  in  the  telluric  field.  The  linear  relation¬ 
ships  of  equation  (1)  will  be  true  for  changes  in 
10 

the  field.  Thus,  we  may  write,  expressing  the 
differences  in  the  field  quantities  as  increments: 


Z^  X  =  a  +  b  ^  y 

(2) 

Z^  Y  =  c  ^x  +  d^y 

Equations  (2)  are  similar  to  equations  (1)  and  as 
before,  A,  the  determinant  of  the  coefficients,  is 
unique  at  a  given  location.  The  elements  of  the 
determinant  are  obtained  by  solving  equations  (2), 
using  as  values  of  the  increments,  changes  in  the 
field  quantities  measured  during  some  arbitrary 
time  interval. 

Each  set  of  corresponding  base  and  field 
station  readings  must  be  standardized  so  that  the 
data  from  various  stations  may  be  compared.  There 
will  be  two  sets  of  readings  corresponding  to  the 
quantities  in  equation  (2):  the  readings  at  the 
base  station  and  the  readings  at  the  field  station. 
To  standardize  the  two  sets  of  readings,  it  is  most 


10 


Ibid. .  p.  52 
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convenient  to  normalize  the  base  vectors  so  that 
when  represented  vectorially,  the  end  ^points  lie  on 
a  unit  circle.  To  do  this,  we  set 


2  7  9 

A  z  =  +  Ay^ 

Then  let 

X1  =  ^X 

-A  Z 

Y  1  =  ^  Y 

a  z 


(3) 


x1  =  ^  x 
Az 


y 1  =  a  y 

a  z 

In  terms  of  these  new  quantities,  equation 
(2)  becomes 


X*  =  ax1  +  by* 
Y1  =  ex1  +  dy1 


(4) 


and 

X*2  +  y‘2  =  1  (5) 

Equation  (5)  indicates  that  the  base  station 
electric  vector  lies  upon  a  unit  circle.  Under  the 
same  operations,  the  field  station  vector  will  have 
a  characteristic  form.  This  can  be  found  by  solving 
equation  (4)  for  x*  and  y\ 


1 


x  *  =  dX1  -  bY1 
ad  -  be 
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(6) 


y*  =  aY1  -  cX* 
ad  -  be 

Combining  equations  (5)  and  (6)  gives 

2  2 

(dX1  -  bY*)  +  (aY1  -  eX*)  =  ad  -  be 
This  is  the  equation  of  an  ellipse.  The 
normal  vector  at  the  field  station  thus  describes 
an  ellipse  when  the  normal  vector  at  the  base  sta¬ 
tion  describes  a  unit  circle «  The  area  of  the 
ellipse  is  S  =  tt  (ad  -  be)  =  7f A,  where  A  is  the 
determinant  of  the  linear  vector  transformation 
describing  the  relation  between  the  telluric  field 
at  the  field  station  and  at  the  base  station.  A 
is  thus  the  ratio  of  the  area  of  the  ellipse  at  the 
field  station  to  the  area  of  the  unit  circle  at  the 
base  station.  Since  the  determinant  A  is  independent 
of  the  choice  of  axes,  the  ellipse  area  is  inde¬ 
pendent  of  the  axes  used  to  measure  the  electric 
field*  Therefore,  the  area  and  orientation  of  the 
field  ellipse  at  any  point  is  independent  of  the 
direction  of  the  electric  field* 

The  form  of  the  field  vector  ellipse  gives 
information  about  the  anisotropic  nature  of  the 
subsurface  electrical  resistance.  In  accordance 
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with  Ohm's  law,  the  greatest  potential  gradient  on 
the  earth's  surface  is  in  the  direction  in  which  the 
current  encounters  the  greatest  resistance.  This 
direction  is  -that  of  the  major  axis  of  the  ellipse* 
Perpendicular  to  this,  the  current  encounters  the 
least  resistance,  and  this  is  the  direction  of  the 
minor  axis  of  the  ellipse*  An  especially  high 
apparent  resistance  would  produce  a  large  ellipse 
area*  Thus,  the  absolute  magnitude  of  the  electric 
vector  may  be  interpreted  as  a  relative  quantity 
expressing  the  apparent  resistance  of  the  soil. 

The  special  structural  case  of  this  research 
is  that  of  two  layers,  a  sedimentary  layer  over¬ 
lying  a  highly  resistant  crystalline  basement.  With 
this  type  of  structure,  the  area  of  the  ellipse  is 
inversely  proportional  to  the  thickness  of  the 
sedimentary  layer.  Thus,  a  profile  or  a  map  showing 
variations  in  the  area  of  the  ellipses  over  some 
area  should  give  a  qualitative  picture  of  the  upper 
surface  of  the  basement  below  that  area. 


16 


c)  Magneto-telluric  Theory 

The  general  purpose  of  the  magneto-telluric 
prospecting  which  was  investigated  is  to  obtain 
structural  information  by  matching  measured  values 
of  the  magnetic  and  telluric  fields  at  the  surface 
with  theoretical  values*  Before  this  can  be  done, 
the  theoretical  relationships  between  the  compon¬ 
ents  of  the  electromagnetic  field  at  the  surface 
must  be  developed*  To  do  this,  certain  properties 
of  both  the  electromagnetic  field  and  the  subsurface 
structure  are  assumed  and  under  these  assumptions 
the  theoretical  relations  are  obtained* 

First,  it  is  assumed  that  the  variations  in 

the  magnetic  field  induce  the  telluric  currents* 

It  is  further  assumed  that  the  induced  current 

sheet  is  uniform  and  in  one  direction  only.  The 

assumption  of  uniformity  is  necessary  and  is 

11 

verified  by  experiment*  The  assumption  that  the 
currents  flow  in  one  direction  is  made  for  conven¬ 
ience  in  development  of  the  theory  and  does  not 
result  in  any  loss  in  generality. 

With  these  assumptions,  the  currents  at  the 

~ 

Schlumberger  and  f^unetz ,  Trans  Am  Geoph 
Union,  3:271  (1939). 
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surface  may  be  expressed  as: 


I 


x 


COS  00  t 


I. 


y 


I 


z 


0 


Then  at  depth  z: 


is  the  conductivity  of  the  soil* 

If  p  (depth  of  penetration)  is  the  depth  at 
which  the  amplitude  is  reduced  to  of  its  value 
at  the  surface,  then 


1 


P 


YxtptZo 


It  can  be  seen  that  as  to,  the  angular  fre¬ 
quency  increases,  the  depth  of  penetration  becomes 
less.  That  is,  phenomena  of  a  short  period  have 
less  effect  at  depth  than  those  of  a  longer  period* 
This  characteristic  of  electromagnetic  quantities 
is  known  as  the  ’’skin  effect,"  It  has  a  marked 
effect  on  magneto-telluric  relationships,  since 
the  subsurface  effects  of  an  electromagnetic 
phenomenon  decrease  inversely  as  the  period  of  the 

12 

The  theoretical  development  follows  the 
outline  of  Cagniard,  Geophysics .  18:605  (1953) 
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phenomenon.  The  depth  of  electromagnetic  pene¬ 
tration  as  a  function  of  period  is  shown  in  Figure 

1. 

The  relation  between  the  electric  and  magnetic 
fields  at  the  surface  of  the  earth  may  be  developed 
using  the  Hertzian  vector,  II,  which  satisfies 
Maxwell's  equations.  In  terms  of  the  Hertzian 
vector,  Maxwell's  equations  give  for  a  horizon¬ 
tally  stratified  earth: 

V2II  +  =  0  (9) 


H  =  4^0" curl  II 

or  H  =  4  TT  CT  Dll 
c)  z 

Hx  =  Hz  =  0 


(10) 


so  that 


E  =  grad  divll  - 
Ex  =  iix 


Ey  =  Ez  =  0 


(11) 


Since  Ex  is  proportional  to  IIX,  Ex  may  be 
chosen  as  the  Hertzian  vector  without  loss  of 
generality.  Then  (9)  is: 


+  477<T^EX  =  0 


H 


y 


—  o  c)  Bx 
a)  3>z 


and 
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ELECTROMAGNETIC  PENETRATION 


FIGURE  1 

DEPTH  OF  ELECTROMAGNETIC  PENETRATION  AS  A  FUNCTION 
OF  PERIOD  (After  Cagniard ) 
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Ev  has  the  form 

A 


where 


Ex  =  A  exp(a'V<rz)  +  B  exp(-aVFz) 
!  a  =  2  if exp(-  i  -7^-) 


If  the  subsurface  is  homogeneous,  the  relation 
between  the  electric  and  magnetic  vectors  will  be: 


(12) 


In  equation  (12),  the  value  of  the  resistivity, 
f>  ,  is  the  true  resistivity  of  the  subsurface  if 
the  subsurface  is  homogeneous.  For  the  more  general 
case  of  a  stratified  subsurface,  however,  the 
resistivity  of  the  soil  will  increase  with  depth  so 
that  there  will  not  be  a  unique  resistivity.  The 
value  of  p  obtained  from  equation  (12)  is  known  as 
the  "apparent  resistivity." 

Equation  (12)  is  shown  graphically  in  Figure  2. 
In  this  figure,  the  intensity  in  gammas  (one  gamma 
equals  10  ^  oersteds)  corresponding  to  an  electric 
field  of  one  millivolt  per  kilometer  is  shown  for 
various  resistivities. 

In  the  general  case,  the  subsurface  will  be 
horizontally  stratified  with  layers  whose  resistivity 
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- 


c 


I  , 


. 


FIGURE  2 

MAGNETIC  INTENSITY  CORRESPONDING  TO  AN 
ELECTRIC  FIELD  OF  ONE  MILLIVOLT  PER 
KILOMETER  (After  Cagniard) 
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will  increase  with  depth.  Because  of  the  skin 
effect,  the  longer  period  variations  will  he  more 
affected  by  the  higher  resistivity  of  the  deeper 
layers  than  will  the  shorter  period  variations. 

Thus  the  apparent  resistivity  computed  from  varia¬ 
tions  with  long  periods  will  be  greater  than  for 
those  with  short  periods.  That  is,  the  apparent 
resistivity  is  a  function  increasing  with  period 
if  the  underground  is  horizontally  stratified. 

The  particular  case  to  be  dealt  with  here  is 
that  of  two  layers :  a  sedimentary  layer  of  low 
resistivity  above  a  highly  resistive  crystalline 
layer.  This  roughly  corresponds  to  the  condition 
across  the  Canadian  Prairies,  where  several  thou¬ 
sand  feet  of  sedimentary  layers  overlie  the  crystal¬ 
line  Precambrian  basement.  While  it  is  true  that 
the  sedimentary  layers  do  not  have  homogeneous 
electrical  properties,  the  resistivity  of  the 
Precambrian  basement  is  greater  than  the  resistivity 

of  the  overlying  sediments  by  a  factor  of  the  order 

6 

of  10  .  The  greatest  variation  of  resistivity  in 
the  sediments,  on  the  other  hand,  is  usually  less 
than  a  factor  of  102  between  the  upper  and  lower 
sedimentary  layers. 


23 

The  electromagnetic  relations  at  the  surface 


corresponding  to  the  two-layer  case  may  be  deter¬ 
mined  by  using  the  values  of  (13)  and  (14)  and 
solving  for  two  layers.  In  the  upper  layer: 

Ex  =  A  exp(aY^Tz)  +  B  exp(-aVoT  z) 

U5) 

Hy  =  exp(^~V2^VT)  (-A  exp(aY^z)  +  B  exp  (-a^^Tz ) 

In  the  lower  layer: 

Ev  =  exp  (-a^vz) 

x  , _  (16) 

Hy  =  exp  ( ^  2  ^  T )  exp  (-al^z) 

The  subscripts  1  and  2  signify  the  upper  and 
lower  layers  respectively.  The  conditions  of 
continuity  across  the  boundary  enable  the  factors 
A  and  B  to  be  determined: 

A  ~~  exp  ( -ah  ( )) 

A  =  2^ 

B=  •+  exp  -  <IT%) 

The  fields  at  the  surface  can  be  expressed  in 
terms  of  the  depth  of  penetration,  p,  for  layer  1 
by: 

Ex  =  M  exp-  i  <j>  ( 18  ) 

Hy  =1H7t  N  expi.(2r-^  )  (19) 
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where 

M  cos  <f>  =  (“cosh-^  + j^sinh  ^-)cos  -jjj- 
M  sin?  =  (-^  sinh-~  +  ^cosh  )  sin  A 
N  cos  vf  =  (-^  sinh-p  +  p;cosh ^)cos  ~ 

N  sin'/'  =  (-^  cosh-^+^.sinh-^)sin  ^ 

The  ratio  of  the  electric  component  to  the 
magnetic  component  is: 

Jx_  =  - — — -  JL  exp -i(£  +  4-\i>)  (21) 

Hy  T  N 

The  phase  difference  between  the  telluric  and 
the  magnetic  components  is  given  by  equation  (21) 
and  is  (X-f^-^)*  The  apparent  resistivity,  A,, 
is : 

1  a  JN  1 

Y2C^  T  N  iTr,  T 

or, 


Equations  (20)  and  (22)  may  be  used  to  plot 
theoretical  curves  of  apparent  resistivity  against 
time  for  given  values  of  the  resistivity  of  the  two 
layers  and  the  thickness  of  the  upper  layer » 
Equation  (21)  gives  a  similar  theoretical  curve  for 
the  phase  difference  between  the  telluric  and 
magnetic  components. 


. 


■ 
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A  possible  method  of  magneto-telluric  pros¬ 
pecting  would  be  to  measure  two  horizontal  telluric 
and  magnetic  components  at  right  angles.  Using  the 
values  of  these  components,  which  would  correspond 
to  Ej^  and  H  ,  the  apparent  resistivity  could  be 
computed  using  equation  (14),  The  values  of  appar¬ 
ent  resistivity  obtained  could  then  be  plotted 
against  T,  the  period  of  the  disturbance  used  in 
the  calculations,  and  the  result  could  be  compared 
with  the  theoretical  curves.  The  details  of  curve 
matching  will  be  given  in  detail  in  the  chapter  on 
Methods  of  Analysis, 


CHAPTER  THREE 


FIELD  METHODS 

a)  Telluric  Prospecting 

In  the  summer  of  1956  >  telluric  measurements 
were  made  at  a  number  of  stations  along  a  profile 
across  the  Canadian  Prairies.  (See  Figure  3)*  The 
stations  were  chosen  to  be  on  the  same  geomagnetic 
latitude  as  the  base  station,  the  Dominion  Observa¬ 
tory  at  Meanook,  Alberta.  At  the  westernmost 
station,  Clairmont,  Alberta,  the  top  of  the  granitic 
Precambrian  basement  is  11,500  feet  below  the 
surface.  This  depth  decreases  towards  the  east  as 
the  thickness  of  the  overlying  sediments  becomes 
less,  and  at  Portage  La  Prairie  ,  Manitoba,  the  top 
of  the  Precambrian  is  2000  feet  below  the  surface. 
The  purpose  of  the  telluric  measurements  was  to 
attempt  to  relate  the  variations  in  the  character 
of  the  telluric  currents  along  the  profile  to  the 
change  in  the  thickness  of  the  sedimentary  layers. 

The  recorders  used  to  measure  the  telluric 
field  were  two  Varian  type  G-10.  Recorders  of  this 
type  work  on  a  null-balance  principle;  that  is,  at 
balance  they  draw  no  current.  The  input  voltage  to 
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LOCATION  OP  TELLURIC  AND  MAGNETO-TELLURIC  FIELD  STATIONS 
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the  recorder  is  compared  with  a  standard  cell 
voltage,  and  the  difference  voltage  is  amplified 
by  a  d-c  amplifier.  The  d-c  amplifier  drives  a 
servo-mechanism  which  moves  a  recording  pen  on  a 
moving  paper  record.  The  response  time  of  the 
recorders  (the  time  required  for  the  recording  pen 
to  travel  the  full  width  of  the  scale)  was  one 
second.  In  order  to  achieve  critical  damping  of 
the  system,  the  input  resistance  must  be  less 
than  0.1  megohm.  Usually,  the  signal  resistance 
between  the  electrodes  is  less  than  this  value, 
but  in  some  locations,  notably  those  which  were 
particularly  dry,  this  resistance  was  exceeded, 
and  usable  records  could  not  be  obtained.  In 
these  instances,  when  the  system  cannot  be  critic¬ 
ally  damped,  the  recording  pen  will  not  follow  the 
signal,  and  produces  a  record  with  square-shaped 
variations  as  shown  in  Figure  4. 

The  sensitivity  of  the  recorders  was  stable 
enough  to  meet  the  requirements  of  field  operations. 
The  standard  cells  in  the  recorder,  mercury  hearing 
aid  type  cells,  have  a  very  low  drop  in  E.  M.  F. 
throughout  their  life.  They  need  to  be  replaced 
about  every  four  months,  and  the  sensitivity  during 
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TELLURIC  RECORD  WITH  INPUT  IMPEDANCE  TOO  HIGH 
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this  period  is  fairly  stable. 

The  recorders  had  two  available  paper  speeds: 
four  inches  per  minute  and  four  inches  per  hour.  It 
was  found  that  while  the  faster  paper  speed  gave  the 
high-frequency  variations  which  the  lower  speed 
obscured,  the  difficulties  in  accurately  timing  the 
records  made  the  fast  paper  speed  impractical. 

Thus,  all  the  records  were  obtained  at  a  paper 
speed  of  four  inches  per  hour*  One  tef  the  recorders 
was  used  for  the  North-South  component,  and  the 
other  for  the  East-West.  Timing  marks  were  pro¬ 
duced  on  the  records  by  applying  a  small  voltage  to 
the  recorders  at  regular  intervals. 

The  records  at  the  Meanook  Observatory  base 
station  were  taken  on  Brown  recorders.  These 
recorders  have  the  advantage  that  the  two  telluric 
components  are  recorded  on  the  same  paper  record, 
which  reduces  the  error  in  timing. 

The  input  terminals  of  the  recorders  were 
connected  to  the  electrodes  which  were  placed  1000 
feet  apart.  The  usual  field  practice  was  to  install 
the  four  electrodes  five  hundred  feet  from  the 
recording  truck,  which  was  at  the  centre  of  the 
cross  formed  by  the  cables  going  out  to  the  elect- 
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rodes  The  cables  used  were  five  hundred  feet 
lengths  of  number  18  rubber  insulated,  flexible 
copper  wire*  These  cables  were  easily  laid  out 
and  wound  in  at  each  location  by  the  hand  held 
reel  shown  in  Figure  5* 

The  electrodes  themselves  were  copper  rods, 
immersed  in  an  electrolyte  of  copper  sulfate,  which 
was  contained  in  unglazed  porous  pots  five  inches 
deep*  In  order  to  maintain  a  uniform  condition  of 
the  electrolyte,  the  solution  was  kept  saturated* 

This  type  of  electrode  eliminates  most  of  the 
contact  potential,  described  in  the  previous  chapter* 
In  practice,  it  was  found  that  there  was  a  smaj.1, 
slightly  varying  contact  potential  as  the  electro¬ 
lyte  in  the  electrodes  seeped  out*  This  was  poss¬ 
ibly  due  partly  to  a  change  in  the  damping  char¬ 
acteristics  of  the  recorders  as  the  resistance  of 
the  electrodes  varied*  Figure  5  shows  an  electrode 
installed  in  the  ground* 

The  instruments  were  operated  on  115  volts 
supplied  by  an  ATR  inverter  operating  from  storage 
batteries*  The  advantage  of  the  inverter-type 
supply  is  that  it  is  cheap.  However,  it  draws  a 
heavy  current  from  the  storage  batteries,  which 
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TYPICAL  ELECTRODE  INSTALLATION 
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makes  the  voltage  output  drop  with  use.  This 
variation  in  line  voltage  affected  the  paper  speed 
and  sensitivity  of  the  recorders.  Another  dis¬ 
advantage  of  the  inverter-type  supply  is  that  it  is 
not  recommended  for  long  period  operation.  Thus, 
the  recordings  with  the  inverter  power  supply  were 
not  continued  for  more  than  five  hours. 

Figure  6  shows  the  recording  apparatus  in¬ 
stalled  in  the  recording  truck  during  a  recording 
period.  A  voltmeter  was  used  to  make  periodic 
checks  on  the  line  voltage,  which  was  adjusted  by 
a  control  on  the  inverter.  The  short-wave  radio 
was  used  to  receive  time  signals.  The  records  were 
timed  as  accurately  as  possible  to  enable  them  to  be 
compared  with  the  corresponding  records  obtained  at 
the  Meanook  Observatory, 
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35 


b)  Magneto- telluric  Prospecting 

In  the  spring  of  1957,  magneto-telluric 
recordings  were  made  at  five  locations  in  Alberta 
and  Saskatchewan,  shown  in  Figure  3*  These  stations 
were  chosen  to  give  a  variation  in  the  depth  of  the 
resistive  basement  below  the  surface*  This  depth  of 
basement  varied  from  12,000  feet  at  McLeod  Valley, 
Alberta,  to  1,600  feet  at  Molanosa,  Saskatchewan* 

In  the  design  of  apparatus  to  measure  the 
magneto-telluric  components,  it  is  necessary  to 
have  a  magnetic  detector  which  will  record  varia¬ 
tions  small  enough  to  compare  with  the  telluric 
variations.  In  Table  I,  the  magnetic  field  in 
gammas  (one  gamma  =  10“^  oersteds)  corresponding 
to  a  telluric  field  of  one  millivolt  per  kilometer 
is  given*  Under  typical  conditions,  the  resistivity 
of  the  soil  would  be  about  fifty  ohm-meters,  and  the 
telluric  variation  about  ten  millivolts  per  kilo¬ 
meter*  Lrom  Table  I,  it  can  be  seen  that  the 
magnetic  variation  corresponding  will  have  an 
amplitude  of  eleven  gammas  for  a  variation  with  a 
period  of  five  minutes.  Thus,  an  adequate  recording 
magnetometer  should  be  able  to  give  a  sensitivity 


(After  Gagniard) 
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which  can  be  read  at  least  to  the  nearest  five 
gammas . 

In  order  to  measure  variations  of  the  order 
of  at  least  five  gammas,  we  must  be  able  to  measure 
variations  which  are  only  a  tiny  fraction  of  the 
total  magnetic  intensity.  The  horizontal  compon¬ 
ent  of  the  magnetic  intensity  at  Meanook,  for 
instance,  is  13,000  gammas,  so  that  five  gammas  is 
one  part  in  twenty-five  hundred.  A  small  variation 
of  this  magnitude  may  be  measured  with  a  flux-gate 
magnetometer ,  which  cancels  out  all  but  a  small 

portion  of  the  total  magnetic  field  by  means  of  a 

13 

steady  counter  magnetic  field.  The  variations 

are  then  large  enough  Relative  to  the  remaining 

portion  of  the  magnetic  field  to  be  measured. 

The  flux-gate  magnetometer  system  consists 

of  a  detecting  head,  a  phase  sensitive  detector 

and  amplifier,  a  bucking  circuit  and  a  calibrating 

circuit.  The  phase  sensitive  detector  circuit  is 

shown  in  Figure  ?•  It  is  adapted  from  one  described 

1  A 

by  Meek  and  Hector. 

Wyckoff,  Geophysics .  13:182  (1948). 

14  Meek  and  Hector,  Gan.  Jour,  of  Physics. 
33:364  (1955) 
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The  detecting  head  of  the  magnetometer,  shown 
in  Figure  8,  consists  of  two  matched,  mu-metal 
cores  around  which  have  been  wound  equal  and 
oppositely-acting  primary  windings,  A  single  sec¬ 
ondary  winding  is  wound  outside  of  the  primary 
coils,  which  are  in  series.  The  push-pull  oscill¬ 
ator,  shown  in  Figure  7,  produces  a  signal  of 
about  five  hundred  cycles  per  second  which  is 
applied  to  the  primary  coils  of  the  detecting  head. 
When  no  magnetic  field  is  present,  no  signal  is 
detected  by  the  secondary,  as  the  primary  coils  are 
similar  and  oppositely  wound.  When  an  external 
field  does  exist,  the  primary  applied  magnetic 
field  is  shifted  with  respect  to  the  hysteresis 
loop  of  the  mu-metal  cores.  As  the  cores  operate 
in  opposite  phases,  the  resultant  magnetic  flux 
produces  a  signal  which  is  detected  in  the  second¬ 
ary  at  double  the  primary  frequency.  This  signal 
is  amplified  and  compared  with  double  the  oscill¬ 
ator  frequency  in  the  phase  sensitive  detector. 

The  result  is  amplified  in  the  d-c  amplifier  and 
used  to  operate  a  one  milliampere  recording  meter. 

The  detecting  head  is  mounted  in  the  centre 
of  a  pair  of  Helmholtz  coils,  which  are  used  to 


>  t 

, 


• 

. 

' 

, 


40 


DETECTING  UNIT,  FLUX-GATE  MAGNETOMETER 


calibrate  the  magnetometer.  The  helmholtz  coils 
can  be  seen  in  Figure  8,  A  small  amount  of  current 
was  applied  to  the  calibrating  helmholtz  coils 
about  every  half  hour  during  operation.  Since  the 
field  in  the  region  of  the  detecting  head  was 
known  for  a  given  current,  the  record  could  be 
calibrated.  The  calibrating  mark  also  served  as  a 
time  mark. 

A  further  winding  around  the  detecting  head 
applied  a  bucking  current  which  was  adjusted  so 
that  only  a  small  portion  of  the  earth1 s  field 
affected  the  magnetometer.  Only  a  few  milliamperes 
were  necessary  for  this  purpose.  The  bucking 
current  was  adjusted  from  time  to  time  so  that  the 
recording  magnetometer  would  not  go  off  scale  with 
large  variations  in  the  magnetic  field.  The  bucking 
current  was  supplied  by  storage  batteries  and  did 
not  vary  noticeably  during  the  four  or  five  days 
recording  at  each  station* 

The  d-c  output  of  the  unit  is  designed  to  be 
used  with  a  one  milliampere  recording  meter.  By 
using  a  voltage  dividing  system,  the  output  could 
also  be  used  with  a  Varian  G-10  recording  potentio¬ 
meter  as  was  used  for  the  telluric  measurements. 
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It  was  found  that  the  Varian  recorder  gavd  a  better 
record  since  the  high  frequency  response  of  the 
unit  reproduced  high  frequency  disturbances  well* 
The  recording  milliammeter ,  which  draws  a  current, 
does  not  follow  high  frequency  variations  accur¬ 
ately*  Most  of  the  magnetometer  records  were 
obtained  from  the  recording  milliammeter,  due  to  a 
failure  of  the  Varian  recorder# 

The  magneto-telluric  recording  apparatus  was 
installed  in  a  panel  truck,  and  was  powered  by  a 
fifteen  hundred  watt,  two-cycle  gasoline  generator. 
The  output  from  the  generator  was  not  steady  enough 
to  operate  the  recording  apparatus,  even  though 
a  regulating  system  was  built  into  the  output  of 
the  generator*  In  order  to  give  more  stability  to 
the  line  voltage,  the  output  of  the  generator  was 
applied  to  a  constant  voltage  transformer.  In 
addition,  the  magnetometer  circuit  was  provided 
with  an • electronic  voltage  control  circuit#  The 
high  voltage  for  the  magnetometer  was  supplied 
by  a  General  Radio  Type  12 03 A  Unit  Power  Supply# 

In  operation,  the  detecting  unit  of  the 
magnetometer  was  set  up  about  two  hundred  feet  away 
from  the  recording  truck.  Electrodes  to  measure 
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the  telluric  field  were  installed  at  right  angles 
to  the  direction  of  the  magnetometer  and  five 
hundred  feet  on  each  side.  That  is,  the  magneto¬ 
meter  detector  was  placed  at  the  centre  of  the 
telluric  cables  and  at  right  angles  to  the  direction 
of  the  cables. 

The  component  of  the  telluric  field  was 
recorded  in  the  same  manner  as  described  for  the 
telluric  prospecting.  However,  the  sensitivity  of 
the  Varian  recorder  was  increased  to  ten  milli¬ 
volts  full  scale,  which  provided  a  more  usable 
record  than  the  previous  sensitivity  of  fifty 
millivolts  full  scale.  Since  this  higher  sensitiv¬ 
ity  was  too  great  to  record  on-scale  during  periods 
of  intense  electomagnetic  activity,  a  two-hundred 
ohm  attenuator  was  used  to  divide  the  signal  input 
by  a  convenient  factor.  Using  this  attenuator,  the 
source  resistance  could  be  standardized  to  two 
hundred  ohms  so  that  the  recorder  could  always  be 
critically  damped. 


CHAPTER  POUR 


Methods  Of  Analysis 

a)  Analysis  of  Telluric  Records 

The  area  of  the  ellipse  at  a  field  station 
is  obtained  from  a  record  such  as  the  one  shown  in 
Figure  9*  Only  a  single  component  is  shown  in  this 
figure;  the  corresponding  record  of  both  components 
at  the  base  station  is  shown  in  Figure  10.  As  was 
mentioned  in  the  theory,  changes  in  the  telluric 
vector  are  used  in  the  analyses  of  the  records. 

The  changes  during  any  time  interval  are  given  a 
sign  depending  on  whether  the  telluric  current  is 
increasing  or  decreasing  during  the  interval.  A 
positive  difference,  for  instance,  corresponds  to 
an  increasing  current. 

The  first  step  in  the  analysis  of  the  records 
is  to  read  off  the  changes  in  the  two  components  for 
a  large  number  of  time  intervals.  That  is,  for  each 
time  interval,  &tj_,  the  changes  and  AYi  will 

be  obtained  from  the  field  records  and  the  changes 
and  A  yj_  from  the  base  records.  We  have,  from 
equation  (3): 

2  2  2 
=  ^xi  *  Ayi 
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TYPICAL  TELLURIC  RECORD  OBTAINED  AT  FIELD  STATION 


FIGURE  10 

TYPICAL  TELLURIC  RECORD  OBTAINED  AT  BASE  STATION 
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Then,  as  was  shown  in  the  theory,  the  vectors 
Xj_ 1 ,  and  yi 1  at  the  base  station  form  the  vector  vi  1 
which  lies  on  a  unit  circle.  The  vectors  1  and 
Y^  1  form  the  vector  VA  which  lies  on  an  ellipse  at 
the  field  station. 

In  order  to  obtain  a  numerical  method  for 
the  computation  of  the  ellipse  areas,  let  us  consider 
the  difference  vectors  for  two  different  intervals, 
t^  and  t2*  ^ 

Then,  at  the  base  station: 

=  vx 

+  A  yv>  =  v2 

And,  at  the  field  station: 


15 

The  method  of  analysis  follows 
Porstendorf er ,  Tellur Ik.  pp.  61-67,  Berlin,  195^- 


Now, 


According  to  equation  (1)  then, 


=  a  ^x^ 

^Y1 

=  c^x^ 

= 

2 

=  c^T2 

^2 


Forming  the  vector  product  of  the  two  vectors 
x  V2I  =  -  &X2  Al! 

=  (ad  -  be)  (.^  x^y-g  -  Ax2^y1) 

S 1  =  S  »  ad  -  be 

IT 


^Xi/^Y2  -  ax2ay1 

axi^y2  „  ^X2 ayi 


(21) 


Therefore,  in  terms  of  the  vector  products  of  the 
base  and  field  difference  vectors,  we  can  express 
the  ratio  of  the  ellipse  area  of  the  unit  circle  as 
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Using  this  equation,  the  area  of  the  ellipse 
may  be  determined  except  for  the  factor  TT  by 
forming  the  quotient  of  the  vector  product  of  the 
two  difference  vectors  at  the  field  and  base  stations. 

The  vectors  v^  and  may  be  plotted  to  give 
a  graphical  picture  of  the  relation  between  the  unit 
circle  and  the  field  ellipse.  It  is  difficult  to 
determine  the  area  of  the  ellipse  by  this  method 
unless  a  large  number  of  difference  vectors  are 
analyzed. 

The  computations  involved  in  determining  the 
ratio  of  the  ellipse  area  to  the  unit  circle  are 
quite  lengthy  if  the  ratio  is  to  be  determined 
reliably.  A  simpler  method  of  analysis  is  to 
compute  the  ratios  of  the  absolute  values  of  the 
corresponding  difference  vectors  at  the  base  and 
the  field  stations.  That  is,  the  ratios 
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I  vi| 

|Vi\ 


(22) 


may  be  determined  for  a  large  number  of  time  inter¬ 
vals,  and  an  average  taken.  The  average  value  is 
a  measure  of  the  apparent  specific  resistivity  of 


the  soil 
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b)  Analysis  of  Magnet o-telluric  Records 

The  equations  developed  in  the  magneto- 
telluric  theory  are  based  on  the  assumption  that 
the  electromagnetic  variations  are  composed  of  a 
spectrum  of  harmonic  terms*  In  order  to  compare 
the  magnetic  variations  of  a  specific  frequency, 
therefore,  it  is  first  necessary  to  separate  the 
harmonic  components  of  the  record*  An  harmonic 
analysis,  or  Fourier  analysis,  enables  a  cyclic 
function  to  be  expressed  as  a  sum  of  sine  terms 
with  respective  phase  angles.  Then  the  amplitude 
ratios  of  corresponding  harmonics  of  the  telluric 
and  magnetic  records  might  be  used  to  give  the 
apparent  resistivity  in  accordance  with  equation 
(12)*  In  general,  it  is  not  possible  to  separate 
the  corresponding  component  harmonics  by  a  simple 
inspection  of  the  record,  since  two  series  composed 
of  harmonic  terms  out  of  phase  with  each  other 
will  not  have  readily  identifiable  similarities. 

On  the  records  obtained,  it  was  noticed  that 
there  were  some  portions  which  were  cyclic.  That 
is,  similar  features  were  repeated  on  the  record  at 
regular  intervals.  These  particular  portions  were 
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analyzed  according  to  standard  methods  described 

16 

by  Conrad  and  Pollack,  to  give  the  harmonic 
components . 

The  cycles  on  the  record  chosen  to  be  analyz¬ 
ed  are  divided  into  equally  spaced  time  intervals 
(usually  minutes  or  half-minutes . )  The  value  of 
the  ordinate  at  each  interval  is  read  from  the 
records.  These  values  are  taken  over  as  long  an 
interval  of  the  records  as  is  possible,  providing 
that  the  record  is  cyclic  over  the  entire  interval. 
The  more  cycles  over  which  the  data  is  read,  the 
less  is  the  effect  of  random  fluctuations. 

The  data  read  from  the  records  is  entered  on 

a  “folded"  form.  “Folding"  the  data  will  reduce  it 

n  17 

to  one  average  cycle.  Suppose  that  observations 
have  been  made  for  n  cycles.  Each  cycle  will  have 
the  same  number,  r,  of  individual  observations  so 
that  altogether  there  will  be  N  =  r  x  n  values  of 
the  ordinates.  We  arrange  the  N  data, 

yGi  y±f  y2>  •♦••••••  •  yN-1 

in  a  table  with  r  columns.  Then  the  original  data 


Conrad  and  Pollack,  Methods  in  Climatology. 
p.119,  Cambridge,  1950. 

17  Ibid.,  pp.  355-58. 
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y  ^  appear  folded  into  n  rows  of  r  values  each.  In 
the  last  row,  the  averages,  Y^,  for  each  column  are 
entered.  This  last  row  of  averages  is  the  row  of 
the  ordinates  of  an  average  cycle  of  the  distur¬ 
bance  chosen.  The  average  row  preserves  the  basic 
period  and  all  its  submultiples  and  weakens  the 
random,  low-persistance  fluctuations. 

If  the  average  cycle  thus  obtained  were 
periodic,  the  ordinates  at  each  end  of  the  cycle 
would  be  the  same.  Usually,  though,  this  is  not 
the  case,  because  of  the  effect  of  non-cyclic  or 
long  term  changes.  Before  the  data  is  suitable  for 
an  harmonic  analysis,  this  long  term  change  must  be 
removed  from  the  cycle.  This  is  done  by  applying  a 
small  linear  correction  to  each  value  of  the  ord¬ 
inates  so  that  the  two  end  points  are  the  same. 

The  center  ordinate  is  used  as  the  mid-point  of  the 
correction,  and  is  left  unchanged. 

The  non-cyclic  correction  makes  the  assumption 
that  the  drift  is  linear  during  the  interval  of  the 
cycles  analyzed.  This  would,  in  general,  not  be 
true.  However,  the  non-linear  error  is  small  if 
several  cycles  are  analyzed.  The  linear  correction 
serves  at  least  to  eliminate  the  major  portion  of 
the  change  which  would  otherwise  produce  erroneous 
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results  in  the  harmonic  analysis* 

The  method  of  harmonic  analysis  employed  uses 
twelve  ordinates  and  gives  five  harmonics  with 
their  respective  phase  angles*  The  analysis  is 
most  conveniently  carried  out  on  a  prepared  form, 
such  as  the  one  adapted  from  Whittaker  and  Robin¬ 
son  and  shown  on  Page  54*  ^ 

The  twelve  equally-spaced  ordinates  are 
taken  from  the  average  cycle.  In  the  cases  where 
the  number  of  ordinates  for  which  the  cycle  is 
computed  is  not  a  multiple  of  twelve,  the  average 
cycle  was  plotted  and  twelve  ordinates  were  read 
off  its  graph. 

Each  pair  of  corresponding  magnetic  and 
telluric  records  were  treated  in  a  parallel  manner. 
That  is,  the  data  was  read  from  both  records  for  the 
same  periods,  and  the  same  ordinate  spacing  was 
used.  The  results  of  the  analysis  give  the  amp¬ 
litude  and  phase  angles  of  the  first  five  harmonics 
of  both  fields.  The  amplitude  ratios,  (i.e.,  Ex) 


for  each  of  the  harmonics  is  then  readily  determined. 


Whittaker  and  Robinson,  The  Calculus  of 
Observations .  pp.  267-71,  London,  1944. 
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Using  equation( 12 ) ,  the  amplitude  ratio  will  give  a 
value  of  apparent  resistivity  for  the  location. 

This  value  of  apparent  resistivity  will  vary  with 
the  period  of  the  harmonic  used  if  the  subsurface 
is  not  homogeneous. 

As  was  mentioned,  equations  (18),  (19),  and 

(20)  may  be  used  to  construct  theoretical  curves, 

19 

such  as  given  by  Gagniard,  Two  sets  of  curves 
may  be  plotted  from  the  data:  1)  the  apparent 
resistivity  as  a  function  of  the  resistivity  of  the 
upper  and  lower  layers  and  the  thickness  of  the 
upper  layer  (Figure  11),  and  2)  the  phase  difference 
between  the  electric  and  magnetic  fields  as  a 
function  of  resistivity  of  the  layers  and  the  depth 
of  the  upper  layer  (Figure  12). 

The  theoretical  curves  plotted  in  Figures  11 
and  12  are  computed  for  an  upper  layer  with  resis¬ 
tivity  of  ten  ohm-meters,  and  a  depth,  h,  of  one 
kilometer.  The  set  of  curves  is  computed  for  var¬ 
ious  values  of  the  resistivity  of  the  second  layer, 
Pi  .  To  apply  these  curves  to  a  particular  case 
where  the  resistivity  and  thickness  of  the  upper 
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Cagniard,  Georfeavsics .  18:605-35  (1953)* 
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APPARENT  RESISTIVITY  FOR  TWO-LAYERED  SUBSURFACE 


=  <t>  +•  v  +  rr/4 
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PHASE  DIFFERENCE  BETWEEN  COMPONENTS  FOR  TWO-LAYERED  SUBSURFACE 
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layer  will  not  in  general  be  the  same  as  for  the 

theoretical  curves,  the  ordinate  of  the  theoretical 

apparent  resistivity  curves  (Figure  11)  must  be 
h2 

multiplied  by  .  Since  the  scales  are  logar- 

r «  c. 

ithmic,  both  these  factors  may  be  applied  to  the 
theoretical  curves  by  a  shift  of  the  axes. 

To  use  the  computed  values  of  the  apparent 
resistivity  or  phase  difference,  the  values  are 
plotted  on  transparent  paper  against  period  to  the 
same  scale  as  the  theoretical  curves.  The  calcul¬ 
ated  values  on  the  transparent  paper  are  then  laid 
over  the  theoretical  curves.  The  shift  of  the 
origin  necessary  to  match  the  curves  determines 
the  value  of  the  resistivity  of  the  second  layer 
and  the  depth  of  the  upper  layer. 

The  depth  of  the  upper  layer,  h,  may  be  read 
from  the  intercept  of  the  curve  of  the  observed 
data  and  the  time  axis.  This  time  intercept,  Tc, 
has  a  value : 

Tc  =  — 6^ -id  .  (23) 

10P1 


If  a  value  of  £>{  ,  the  resistivity  of  the 
upper  layer  is  aasumed,  the  depth  of  the  upper 
layer  may  be  determined  from  equation  (23)* 


CHAPTER  FIVE 


ANALYSIS  AND  INTERPRETATION 

a)  Analysis  of  a  Typical  Telluric  Record 

In  this  section,  a  typical  telluric  record 
will  be  treated  to  show  the  method  of  analysis  in 
detail.  The  record  chosen  was  that  obtained  at 
Clairmont,  Alberta.  Eleven  usable  time  intervals 
were  picked  to  give  the  values  of  the  difference 
vectors  from  the  field  record.  The  time  intervals 
are  only  11  usable"  if  the  records  at  both  the  base 
and  field  stations  are  clear  and  distinct.  Due  to 
calibration  marks,  off-scale  surges  of  the  field, 
and  other  disturbances,  there  were  many  portions  of 
the  records  which  could  not  be  used  for  an  analysis. 

The  difference  quantities  corresponding  to 
the  field  changes  for  the  chosen  intervals  are  shown 
in  Table  II.  The  ordinates  read  off  the  records 
were  first  multiplied  by  suitable  factors  to  give 
the  results  in  units  of  millivolts  per  kilometer, 
which  are  the  units  usually  used  in  telluric 
measurements.  The  amplitude  of  the  base  station 
difference  vector  was  then  computed  and  the  differ¬ 
ence  quantities  were  all  divided  by  this  factor  to 
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normalize  theip.  After  normalization,  the  data  was 
suitable  for  presentation  in  graphical  form.  The 
graph  of  the  base  station  unit  circle  and  the  field 
station  ellipse  are  shown  in  Figure  13. 

The  ratio  of  the  ellipse  area  to  the  area  of 
the  unit  circle  was  computed  using  equation  (21). 
That  is, 


z.X1^Y2  -  ^X2  ^Y1 


^x1^y2  -^x2z^y1 


The  data  obtained  from  the  field  and  base 
records,  and  the  results  of  the  computation  of  the 
field  vector  ellipse  areas  are  shown  on  Table  I* 

The  ellipse  area  for  Clairmont,  Alberta  was 
found  to  be  S'  =  0.44,  with  a  standard  deviation 
of  0.26. 

In  addition  to  the  ellipse  areas,  the  value 
of  the  relative  magnitude  of  the  difference  vectors 
was  calculated  for  each  of  the  intervals.  The 
amplitudes  of  the  two  difference  vectors  were 
computed  as: 


Then  the  ratio  of  the  two  vectors, 


S 
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The  amplitude  ratios  for  Glairmont  had  an 
average  value  of  2.57  and  a  standard  deviation  of 
0.75. 


ANALYSIS  OF  TELLURIC  RECORDS,  CLAIRMONT,  ALBERTA 
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b)  Interpretation  of  the  Telluric  Analysis 

Table  III  showsi  the  results  of  the  telluric 
analysis  for  each  station  computed  in  the  manner 
presented  in  part  (a)  of  this  chapter.  These  results 
are  shown  graphically  in  Figures  14  and  15*  Figure 
14  shows  the  variation  of  the  ellipse  areas  along 
the  profile  and  the  corresponding  variations  in  the 
depth  of  basement.  Figure  15  shows  the  variation 
in  the  relative  amplitude  of  the  telluric  field 
along  the  profile. 

The  correlation  coefficients  between  the 

basement  depth  and  the  values  of  the  ellipse  areas 

and  relative  amplitudes  respectively  were  calculated 

according  to  methods  described  by  Conrad  and  Poll** 

20 

gck,  A  correlation  coefficient  of  zero  signifies 
no  correlation  between  two  series,  and  a  coefficient 
of  one  means  that  the  series  are  identical  or  differ 
by  a  constant  factor.  A  coefficient  of  0,5  or  less, 
when  computed  from  only  a  few  values,  has  little 
significance,  Conrad  and  Pollack  suggest  that  the 
correlation  coefficient  should  be  at  least  six  times 

20 

Conrad  and  Pollack,  Methods  in  Climatology f 
fip.  245-247,  Cambridge,  1950 


* 

• 

■ 

S 

, 

~  .  L 


? • 


•  -  .  , 

* 

- 


STATION  DEPTH  OF  ELLIPSE  STANDARD  RELATIVE  STANDARD 

BASEMENT  AREA  DEVIATION  TELLURIC  DEVIATION 

(feet)  AMPLITUDES 


66 


VO 

O- 

co 

00 

VO 

CO 

ON 

o 

CM 

CM 

o 

VO 

O- 

CO 

tH 

• 

• 

• 

• 

• 

• 

• 

« 

• 

O 

o 

o 

o 

o 

o 

o 

o 

00 

ON 

so 

00 

IN* 

00 

-3- 

o 

VO 

VO 

Os 

CM 

C'N 

\0 

o- 

ON 

-3- 

• 

• 

• 

• 

• 

• 

• 

• 

• 

CM 

CM 

o 

tH 

tH 

t — 1 

T — 1 

o 

o 

so 

00 

so 

so 

CS] 

VO 

o- 

o 

VO 

-3- 

tH 

ON 

CM 

• 

« 

• 

# 

* 

• 

* 

* 

• 

o 

o 

o 

rH 

o 

o 

o 

rH 

VO 

-3- 

CM 

o- 

co 

o- 

-3- 

-3" 

CM 

-3- 

CO 

VO 

VO 

vo 

-3- 

CO 

00 

00 

• 

« 

• 

* 

* 

« 

• 

• 

• 

o 

o 

rH 

o 

o 

o 

o 

rH 

VO 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

VO 

CM 

VO 

o 

CM 

o 

o 

VO 

o 

•s 

•H. 

•V 

** 

m 

•v 

rH 

rH 

10 

SO 

o- 

NO 

so 

VO 

CO 

cvj 

£ 

0 

o5 

-p 

0 

rH 

£ 

•H 

i — 1 

o 

£ 

o 

> 

X 

•H 

•H 

£ 

•H 

0  0 

B 

>5 

o 

> 

rH 

w 

o 

rH 

hD-H 

u 

0 

o 

0 

•H 

m 

-p 

x 

05  t* 

•H 

1 — 1 

Ph 

6 

o 

s 

-P  -H 

o5 

1 — 1 

05 

hO 

Jh 

Jh 

u 

$ 

U  0$ 

rH 

o5 

0 

0 

0 

0 

o 

sm 

o  JL, 

o 

> 

s 

> 

> 

PM  P-< 

67 


UJ 

(/> 


-j  cr 

U  <  <£  in  ^  m  cvj  _  o 


ELLIPSE  AREAS  ACROSS  PRECAMBRIAN  PROFILE 


68 


3IWV8d  VI 

aoviMOd 

NHXNVMd 

N01M8CU 

SS083T 


N0HIM83A 

3TIIA3893A 

MOONV3M 

M3IAA3T1VA 

iNOWbnMO 


-sS 

>  3  j  K 

KtjO 

Sip 


O 

< 


o 

o 

o 


§ 

o 

O' 


O  .  x 
o  t-  £ 
O.  u.  «* 


RELATIVE  TELLURIC  AMPLITUDES  ACROSS  PRECAMBRIAN  PROFILE 


69 


its  probable  error  to  be  considered  significant* 

The  correlation  coefficient  between  the  ellipse 
areas  and  the  basement  was  found  to  be  0,52  -  0*12. 

The  coefficient  is  4.4  times  the  probable  error. 

This  value  has  little  significance  according  to  the 
criteria  mentioned,  so  that  the  results  of  the 
ellipse  area  analysis  across  the  profile  would  not 
seem  to  be  useful  as  an  indication  of  basement  depth* 

The  correlation  coefficient  between  the 
relative  amplitudes  and  the  depth  of  basement  was 
found  to  be  0*68  ^  0.16.  As  before,  the  coeffic¬ 
ient  is  4.4  times  its  probable  error*  Because  of 
the  somewhat  larger  coefficient,  however,  there 
seems  to  be  some  correlation  between  the  relative 
amplitude  of  the  telluric  field  and  the  depth  of 
the  Precambrian  formation,* 

Previous  local  surveys  using  the  ellipse 
areas  have  shown  a  significant  relation  between  the 
ellipse  areas  and  subsurface  structure.  There  are 
several  probable  reasons  why  the  present  investig¬ 
ation  has  not  yielded  more  positive  results  with 
this  method. 

The  actual  records,  themselves,  for  instance, 
were  not  well  suited  to  the  method  of  analysis.  The 
sensitivity  of  the  recorders  used  in  the  telluric 
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exploration  was  fifty  millivolts  full  scale-  This 
was  not  sensitive  enough  to  be  able  to  use  the 
average  variations  in  the  telluric  field  for  an 
analysis,  so  that  only  those  variations  which  were 
somewhat  larger  than  average  could  be  used.  At  most 
stations,  consequently,  only  a  few  intervals  were 
analyzed  with  a  consequent  loss  in  the  significance 
of  the  results.  This  problem  was  later  eliminated 
by  increasing  the  sensitivity  of  the  recorders  to 
ten  millivolts  full  scale. 

Another  source  of  error  in  the  telluric 
analysis  was  the  timing  of  the  records,  While  the 
field  station  records  were  timed  about  every  half 
hour  during  the  four  hour  recording  period,  the 
base  station  records  were  only  timed  about  once 
every  day*  Since  the  paper  drive  at  both  stations 
fluctuated  with  line  voltage,  this  introduced  a 
timing  error  which  could  be  quite  substantial  at 
some  points  on  the  records. 

It  can  be  seen  from  Figure  13  that  the 
variations  in  the  telluric  field  during  the  interval 
of  the  recording  at  Clairmont  all  lie  along  a  single 
narrow  band.  It  was  found  that  during  a  four  hour 
recording  period,  the  direction  of  the  telluric 
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vector  did  not  vary  enough  to  give  a  good  graphical 
indication  of  the  field  vector  ellipse.  This 
difficulty  was  evident  at  the  other  field  stations 
as  well,  so  that  it  seems  necessary  to  record  over 
a  longer  period,  say  twenty-four  hours,  in  order  to 
determine  the  ellipse  graphically. 

The  numerical  computation  of  the  field  vector 
ellipse  has  disadvantages  as  well.  The  quotient  in 
equation  (21)  is  the  difference  between  two  quan¬ 
tities.  If  these  quantities  are  somewhat  similar, 
the  quotient  will  be  very  small,  and  the  probable 
error  extremely  high.  In  many  cases,  the  quotient 
was  nearly  zero,  which  yielded  an  ellipse  area  of 
doubtful  significance. 

The  relative  amplitudes  along  the  profile 
seem  to  have  a  better  correlation  with  basement 
than  the  ellipse  areas.  It  seems  more  practical 
to  use  relative  amplitudes  as  a  method  of  analysis 
than  ellipse  areas.  The  ratio  of  the  amplitudes 
if  far  easier  to  calculate  than  the  ellipse  areas, 
since  the  ratio  of  the  amplitudes  may  be  calculated 
readily  from  the  scale  reading  on  the  records. 

While  the  amplitude  ratios  cannot  give  an  indication 
of  the  anisotropic  nature  of  the  soil,  they  are 
superior  to  the  ellipse  areas  in  that  the  possible 
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error  is  far  less,  and  so  the  values  obtained  are 
probably  more  indicative  of  subsurface  resistivity. 

Similarities  in  the  telluric  records  obtained 
at  the  most  distant  stations  along  the  profile 
indicate  the  uniform  nature  of  the  telluric  currents 
over  a  fairly  large  area  on  the  surface  of  the  earth. 
However,  as  the  distance  between  the  base  station 
and  the  field  station,  increased,  there  was  less  and 
less  similarity  between  the  two  sets  of  records. 

Thus,  it  seems  probable  that  there  is  a  practical 
limit  to  the  area  that  can  be  investigated  with  this 
type  of  survey. 

The  poor  correlation  between  the  results  of 
the  two  analyses  and  the  depth  of  basement  could 
possibly  be  due  to  a  combination  of  two  factors. 

At  the  distant  stations,  that  is  those  at  the  eastern 
end  of  the  profile,  extraneous  fluctuations  in  the 
telluric  field  might  produce  a  record  poorly  suited 
to  analysis.  At  the  western  end  of  the  profile, 
where  the  basement  was  considerably  deeper,  the  mean 
period  of  the  field  changes  used  in  the  analysis  was 
not  great  enough  to  provide  electromagnetic  pene¬ 
tration  to  the  depth  of  the  top  of  the  Precambrian 
formation.  Perhaps  if  an  analysis  were  undertaken 


< 

• 

* 

' ' 

1 

using  field  variations  of  a  period  long  enough  to 
penetrate  to  a  suitable  depth,  more  significant 
results  could  be  obtained • 
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c)  Analysis  of  a  Typical  Magneto-telluric  Record 

The  magneto-telluric  record  chosen  to  he 
analyzed  in  detail  was  obtained  at  Meadow  Lake, 
Saskatchewan.  The  original  records  are  shown  in 
Figures  16  and  17* 

The  values  of  the  ordinates  at  half  minute 
intervals  were  taken  from  the  records  for  the  interval 
from  0050  to  0101  GMT,  June  5#  1957*  This  gave 
twenty-two  values  which  were  folded  on  a  table  on 
the  assumption  that  there  were  two  cycles  of  eleven 
ordinates*  One  average  cycle  was  obtained  from  the 
folded  data,  and  small  linear  corrections  were 
applied  to  the  average  ordinates  of  the  telluric 
and  magnetic  records  respectively*  This  small 
correction  was  the  difference  between  the  end  values 
of  the  average  cycle.  The  final  corrected  average 
cycles  of  the  telluric  and  magnetic  components  are 
shown  in  Figure  18* 

Twelve  equally  spaced  ordinates  were  taken 
from  the  average  cycle  of  both  records.  The  values 
of  these  ordinates  were  entered  on  an  harmonic 
analysis  form.  The  results  of  the  harmonic  analysis 


were : 
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SELECTED  TELLURIC  RECORD,  MEADOW  LAKE 


Fhe  Esterline-Angus  Co.,  Inc.,  indi^apolis.  ind  .  u.s.  a 
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SELECTED  MAGNETIC  BECOED,  MEADOW  LAKE 
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FIGURE  18 

AVERAGE  TELLURIC  AND  MAGNETIC  CYCLES, 
0050-0101,  JUNE  5/57 
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Ex  =  43.8  +  22.4  sin(x  +  1?1°)  +  3.5  sln(2x  +  237°) 

+  1.6  sin(3x  -  78°)  +  0.3  sin(4x  +  38°) 

+  0.8  sin(5x  -  12°) 

H  =69.8  +  6,0  sin(x  +  100°)  +  0.6  sln(2x  +  221°) 

+  0,5  sin(3x  -  51°)  +  0.1  sin(4x  +  18°) 

Each  of  the  above  series  is  expressed  in  scale 
divisions,  so  that  the  numbers  used  must  be  multi¬ 
plied  by  the  scale  factors  of  the  recording  instrument 
to  convert  the  units  to  millivolts  per  kilometer 
and  gammas,  respectively. 

The  apparent  resistivity  which  is  to  be 
determined  is  the  value  of  p  in  equation( 12) , 

Equation  (12)  may  be  put  into  a  form  more  convenient 
for  the  calculation  of  the  apparent  resistivity: 

p*.  =  o.2  t  (fx.  )2 

Hy 

In  this  equation,  T  is  in  seconds,  Ex  is  in 
millivolts  per  kilometer,  and  H  is  in  gammas.  The 
apparent  resistivity  will  be  given  in  ohm-meters, 

The  apparent  resistivity  calculated  from  the  first 
harmonic,  which  has  a  period  of  5*5  minutes,  or 
330  seconds,  will  be: 

Po<.  =  0,2  x  330  x  22,4  x  16,4  ohm-meters 

6.0  X  3.4 


Y 
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Pol  =  2,1  x  10  ohm-meters, 

(The  numbers  "16,4"  and  “3. 4"  in  the  above 
example  are  the  scale  constants  of  the  recorders,) 

The  results  of  the  computations  with  the  other 
harmonics  are: 

Second  Harmonic:  T  =  I65  seconds 

.  .  4 

Pol  =  2,6  x  10  ohm-meters 
Third  Harmonic:  T  =  110  seconds 

Po.  =  5*3  x  10^  ohm-meters 

Fourth  Harmonic:  T  =82.5  seconds 

3 

P^  =  3.5  x  10  ohm-meters 
The  phase  differences  betwTeen  the  harmonics 
are  easily  calculated  by  taking  the  difference  between 
the  corresponding  phase  angles.  Thus,  the  phase 
differences  are: 

First  Harmonic:  71° 

Second  Harmonic:  16° 

Third  Harmonic:  27° 

Fourth  Harmonic:  20° 

These  values  of  apparent  resistivity  and  the 
phase  differences  for  each  harmonic  may  now  be  plotted 
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against  the  period  of  the  harmonic  to  give  observed 
curves  to  be  compared  respectively,  with  the  appar¬ 
ent  resistivity  curve,  Figure  11,  and  the  phase 
difference  curve.  Figure  12.  In  order  to  obtain 
enough  points  for  a  significant  observed  curve, 
several  analyses  of  various  intervals  must  be 
undertaken  to  give  a  large  number  of  points  on  the 


curve 
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d)  Interpretation  of  the  Magneto-telluric  Analysis 

The  results  of  the  apparent  resistivity 
computations  for  Molanosa,  Meadow  Lake,  and  Vega,  are 
shown  in  Figures  20,  21,  and  22  respectively*  It 
Is  evident  that  the  computed  points  do  not  fall  along 
a  curve  similar  to  the  theoretical  curves  shown  in 
Figure  11* 

Although  the  computed  values  are  not  suitable 
for  curve-matching  as  described  in  the  Theory,  there 
seems  to  be  a  general  increase  in  the  apparent 
resistivity  with  period  as  is  predicted  by  the  theory. 

Even  if  the  computed  data  did  lie  approximately 
on  a  curve  similar  to  the  theoretical  curves,  however, 
it  would  not  be  possible  to  determine  the  depth  of 
the  Precambrian  basement  by  curve-matching,  since  the 
range  of  frequencies  used  in  the  analysis  is  very 
small*  There  are  no  values  less  than  one  minute, 
and  only  two  values  greater  than  one  hour.  Conse¬ 
quently,  the  range  for  which  the  computed  curve  is 
determined  is  only  a  portion  of  the  usable  range, 
and  is  not  enough  to  provide  a  unique  match  with  the 
theoretical  curves* 

The  problem  of  finding  a  match  with  the 
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APPARENT  RESISTIVITY,  MEADOW  LAKE 
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theoretical  curves  is  further  exaggerated  by  the 
uncertainty  of  the  resistivity  of  the  upper  layer. 

If  the  resistivity  of  the  upper  layer  were  known 
fairly  accurately,  then  it  would  be  possible  to 
match  even  a  small  portion  of  the  computed  curve 
with  the  theoretical  curve.  If  the  resistivity  of 
the  upper  layer  is  hot  known,  however,  there  is  a 
whole  family  of  possible  theoretical  curves  to 
which  the  computed  portion  may  be  matched. 

The  data  as  presented  is  not  suitable  for 
determining  the  subsurface  properties  as  predicted 
by  Cagniard.  The  limiting  feature  is  the  method  of 
analysis  applied  —  especially  the  harmonic  analysis. 
It  is  not  possible  to  obtain  the  harmonics  over  a 
wide  range  of  periods  or  for  a  large  number  of  cases 
using  the  attempted  method  of  analysis. 

The  small  number  of  points  on  each  of  the 
graphs  of  computed  data  is  due  to  the  difficulty  in 
finding  portions  of  the  record  suitable  to  be  ana¬ 
lyzed,  as  well  as  the  tedious  computations  necessary 
for  the  completion  of  a  single  analysis. 

Although  phase  differences  were  calculated  for 
each  set  of  data,  the  results  were  inconsistent . 

There  was  no  indication  that  the  phase  differences 


. 


, 
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as  computed  could  be  used  for  curve-matching.  The 
computation  of  phase  difference  of  a  short  period 
harmonic  depend  on  an  extremely  accurate  timing  of 
the  records.  It  is  not  believed  that  the  accuracy 
between  the  magnetic  and  telluric  records  was  any 
better  than  plus  or  minus  fifteen  seconds.  This 
magnitude  of  uncertainty  makes  most  phase  difference 
results  doubtful.  For  this  reason,  the  results  of 
the  phase  difference  computations  are  not  shown 
plotted  for  curve-matching. 


CHAPTER  SIX 


CONCLUSIONS 

a)  Summary  of  Findings 

It  was  found  that  the  telluric  prospecting 
as  applied  over  a  large  scale  structural  feature 
showed  little  indication  of  subsurface  structure. 

Two  methods  of  analysis  were  used  with  the 
telluric  data  obtained.  The  first,  the  method  of 
ellipse  areas,  is  favored  by  most  of  the  previous 
investigators  of  telluric  prospecting.  It  was 
found  that  the  method  of  ellipse  areas  has  less 
utility  in  those  instances  where  the  reliability 
of  the  data  is  not  high.  In  these  instances,  an 
analysis  using  the  relative  amplitude  ratios 
provides  a  more  significant  correlation  with 
subsurface  structure  because  of  the  less  complex 
computations  which  are  necessary. 

In  accord  with  the  theoretical  conditions 
21 

predicted  by  Cagniard,  the  telluric  field  mea¬ 
sured  over  a  large  area  showed  a  high  degree  of 
uniformity.  The  condition  of  uniformity  is  a 

21 

Cagniard,  Handbuch  Per  Phvsik.  p.  462, 
Berlin,  1956. 
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necessary  assumption  for  the  development  of  the 
magneto-telluric  theory,  and  also  is  in  agreement 

2 

with  the  results  obtained  by  Sch'lumberger  and  Kunetz. 
There  were  noticeable  fluctuations  present  at  the 
distant  stations  which  wdre  not  present  at  the  base 
station,  however.  As  these  non-uniform  fluctuations 
increase  with  distance  from  the  base,  it  imposes  a 
limit  on  the  scale  of  prospecting  using  the  telluric 
field. 

As  a  general  result,  telluric  prospecting  is 
less  adaptable  to  large  scale  investigations  than[fco 
the  more  local  surveys  where  it  is  possible  to 
obtain  far  greater  control  of  the  records  between 
the  field  station  and  the  base  station. 

It  was  not  possible  to  match  the  computed 
resistivity  curves  with  the  theoretical  curves  in 
accord  with  the  theory  of  Cagniard.  The  major 
limitation  was  that  the  analysis  of  the  magnetic 
and  telluric  records  did  not  cover  a  wide  enough 
range  of  frequencies.  The  computed  data  lay  between 
rather  narrow  frequency  limits  so  that  not  enough  of 

22 

Schlumberger  and  Kunetz,  C..  R.  Acad.  Sci. 
Paris .  223:551  (19^6). 
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the  computed  curve  was  defined  in  order  to  make  a 
unique  match  with  the  theoretical  curves.  For  this 
reason,  the  magneto- telluric  prospecting  investigated 
is  not  practical  unless  a  method  of  analysis  can  be 
devised  which  makes  use  of  a  wider  range  of 
frequencies  • 

The  investigation  of  the  phase-difference 
curve  matching  was  inconclusive  because  the  magneto- 
telluric  records  were  not  timed  accurately  enough 
to  give  reliable  results# 

There  was  no  indication  that  the  theory  as 
proposed  by  Cagniard  is  invalid.  The  restrictions 
of  the  analysis  did  not  permit  enough  of  the  method 
to  be  carried  out  in  order  to  determine  the  limits 
of  Cagniard* s  method  of  subsurface  investigation. 
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b )  Rec  ommendat i ons 

It  is  believed  that  more  positive  results 
might  be  obtained  in  telluric  prospecting  if  the 
timing  of  the  corresponding  records  were  made  more 
accurate*  Also,  it  seems  that  there  is  a  limit  to 
the  scale  of  telluric  prospecting  which  can  be 
undertaken  due  to  the  limits  of  the  uniformity  of 
the  field  over  large  areas. 

It  should  be  made  certain  that  the  variations 
used  in  the  analysis  have  a  period  great  enough  so 
that  their  electromagnetic  penetration  is  deep  enough 
to  be  affected  by  the  sub-surface  structure  inves¬ 
tigated. 


Before  the  magneto -telluric  method  of  pros¬ 
pecting  can  be  made  practical,  a  method  of  analysis 
must  be  devised  which  will  analyse  a  wider  range  of 
frequencies.  It  is  not  believed  that  a  Fourier 
analysis  according  to  standard  methods  will  yield 
a  wide  range  of  frequencies.  Furthermore,  standard 
methods  of  harmonic  analysis  are  only  applicable  to 
phenomena  which  are  periodic.  The  magnetic  and 
telluric  records  obtained  have  few  distinctly  periodic 
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intervals,  and  consist  mainly  of  fluctuating,  low 
persistance  variations. 

The  method  of  analysis  which  seems  to  have  the 
most  promise  of  determining  the  amplitudes  of  the 
component  frequencies  of  the  records  is  that  of  a 
harmonic  analyzer  or  spectrum  analyzer  using 
electronic  methods. 

A  possible  method  of  magneto-telluric  pros¬ 
pecting  would  record  the  corresponding  magnetic  and 
telluric  records  simultaneously  on  magnetic  tape. 
These  could  be  re-played  at  any  desired  speed  in  the 
laboratory.  The  frequency  range  which  may  be 
analyzed  can  extend  over  an  extremely  large  range 
by  choosing  the  speed  at  which  the  records  are  played 
back. 

Unless  the  computed  curves  are  obtained  over 
a  wider  frequency  range,  or  the  resistivity  of  the 
upper  layer  is  known,  the  depth  of  the  lower  layer 
cannot  be  determined.  With  a  wide  frequency  range 
in  the  analysis,  the  resistivity  of  the  upper  layer 
may  be  determined  by  the  high-frequency  variations. 
The  resistivity  of  the  lower  layer  may  be  determined 
by  the  low-frequency  variations,  which  will  pene¬ 
trate  deeper.  If  such  a  wide  range  of  frequencies  is 


* 
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not  available,  it  may  be  difficult  to  investigate 
large  structural  features  by  the  magneto-telluric 
method,  since  it  is  often  difficult  to  ascribe  a 
single  resistivity  which  is  typical  of  the  upper 
sedimentary  layers  as  a  whole. 
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